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ABSTRACT 

Commonplace at every evolutionary stage, Multiple Protostellar Systems 
(MPSs) are thought to be formed through fragmentation, but it is unclear when 
and how. The youngest MPSs, which have not yet undergone much evolution, 
provide important constraints to this question. It is then of interest to disen- 
tangle early stage MPSs. In this letter we present the results of our work on 
VLA1623 using our observations and archival data from the Submillimeter Ar- 
ray (SMA). Our continuum and line observations trace VLAl623's components, 
outflow and envelope, revealing unexpected characteristics. We construct the 
SED for each component using the results of our work and data from literature, 
as well as derive physical parameters from continuum and perform a simple kine- 
matical analysis of the circumstellar material. Our results show VLA1623 to be 
a triple non-coeval system composed of VLA1623A, B & W, with each source 
driving its own outflow and unevenly distributed circumstellar material. From 
the SED, physical parameters and IR emission we conclude that VLA1623A & 
W are Class and I protostars, respectively, and together drive the bulk of the 
observed outflow. Furthermore, we find two surprising results, first the presence 
of a rotating disk-like structure about VLA1623A with indications of pure Keple- 
rian rotation, which, if real, would make it one of the first evidence of Keplerian 
disk structures around Class protostars. Second, we find VLA1623B to be a 
bonafide extremely young protostellar object between the starless core and Class 
stages. 

Subject headings: ISM: individual objects (VLA1623) — ISM: jets and outflows - 
stars: formation — techniques: interferometric 
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Introduction 



Multiple protostellar systems (MPSs) are common among every stage of protostellar 
evolution, but are t ricky to decip her. Fragmentation is believed to be the main mechanism 



for forming MPSs (ITohline 



20021 ). but how and when fragmentation occurs is still a 



major question. To constrain the fragmentation process, information about the system's 



multiplicity, coevality, enviro nment and outf 



suffered dynamical evolution ( jGoodwin et al. 



ows a re needed. Since older systems have 



20071 ). young early stage MPSs, having many 



of the initial conditions almost intact, are the best targets for the task. In this letter we 
present the results of our work on the much debated VLA1623, located in p Ophiuchus and 



generally considered the prototypical Class source flAndre et al 



1993 



, hereafter AWB93). 



19901 ) 



First detected through red- and blueshifted overlapped outflows (I Andre et al. 
VLAl62 3's multiplicity has long been suspected from observati on of multiple continuum 



sources (ILeous et al. 



misalingment ( IDent et al. 



1991 



Strom et al. 



1995 



1995 



Loonev et al. 



Yu & Chernin 



processing jets originating from the same region (IDent et al. 



200(1 hereafter LMW00), outflow 



19971 hereafter YC97) and two discreet 



1995 



Caratti o Garatti et al. 



20061 ). In contrast, the m ultiple continuum sources are argued to co nstitute a monopola r 



jet driven by VLA1623A (IBontemps fc Andre 



Maury et al 



2012 



19971 . hereafter BA97; 



Maury et al. 



2010 



hereafter MOA12). While VLA1623A is agreed to be a real protostellar 



source, and was used by ABW93 as the prime example to establish a new class of YSOs 



'Class 0", there is sti 



(I Ward- Thompson et al. 



1 deb ate on the nature of the other sources, mainly VLA1623B 



2011 



hereafter WKGAll). 



We present the results of our work to determine the mul tiplicity o : 



Ho et al. 



VLA1623 and 



2004 ) observations 



disentangle its outflows. Our own Submillimeter Array (SMA; 
and archival data were used, suplemented with data from literature. Our results show 
VLA1623 to be a triple non-coeval protostellar system and find VLA1623B to be a bona-fide 
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protostar with an extremely young age. 

2. Data 
2.1. Observations 

We observed VLA1623 with SMA at 345GHz in the compact configuration for 8 hours 
on 2007 July 25, with 2GHz correlator mode configured to observe 12 CO (3-2), C 17 (3-2) 
and SO (9-8) simultaneously with 870yum continuum and polarization measurements. All 
windows were set to have 128 channels, providing a velocity resolution of 0.704 km s _1 . 
Quasars 3c273 and 1733-130 were observed as the bandpass and gain calibrators. A 
planet was not observed for unknown operational reasons. This track's original purpose 
was polarization measurements, where lack of flux calibration is not a critical issue. We 
performed standard data reduction for polarization data using MIRIAD. No significant 
polarization detection was made. Flux calibration is done with 1733-130 scaled to 1.31Jy, 
value taken from the SMA Observer Center, with measured uncertainty of 0.035Jy. Imaging 
with natural weighting resulted in a beamsize of 2 / . / 36xl / /41. 

We further observed VLA1623 at 345GHz in the extended configuration with a 3-hr 
filler track on 2012 February 14 using 4GHz correlator mode. Only the continuum emission 
is considered. Calibration was performed on IDL MIR package following the standard 
calibration pipeline. 3c84, 1733-130 and Callisto were used for bandpass, gain and flux 
calibration. Calibrated flux of 1733-130 is 1.0±0.006Jy, in agreement with the flux listed 
on SMA Observer Center. Data imaging and analysis were performed with MIRIAD using 
natural weighting resulting in a beamsize of 1'.'51 xO'/76. 

We compared the two datasets by binning vectorially averaged amplitudes, and when 
the tracks overlap, the flux level was comparable. Uncertainty in flux calibration is typically 
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-10%. 

2.2. Archival Data 

VLA1623 was observed at 230GHz in the very extended configuration on 2007 June 
11 and 2009 July 3. The 2009 data were presented by MOA12. Here we redo the data 
reduction and analysis. The data quality of the 2007 data appears to be better than 
the 2009 data; the system temperatures are 92±24K and 130±39K, and the atmospheric 
opacities are 0.06 and 0.1 for the 2007 and 2009 data, respectively. Both observations used 
the 2GHz correlator mode configured for 12 CO (2-1), C 18 (2-1), 13 CO (2-1) and continuum. 
The three windows with line emission in both datasets were set to have 512 channels each, 
giving a velocity resolution of 0.3kms _1 . The remaining windows were configured to have 
32 channels each. 

For the 2007 data, the quasar 3c454.3 for bandpass, 1514-241 and 1622-297 for gain, and 
Callisto for flux were used in calibration. For the 2009 data, bandpass and gain calibration 
was done with 1924-292 and 1625-254. Flux calibration using Callisto caused the flux of 
VLA1623B to rise 30mJy compared to the 2007 track. We suspect a glitch occured during 
Callisto's observation, so 1625-254 was used as a flux calibrator instead with 0.79Jy taken 
from the SMA Observer Center around the observing date, with measured uncertainty of 
0.025Jy. After scaling, the flux from the two datasets are consistent within 11%. MIRIAD 
was used to separately calibrate each dataset and then combine during imaging assuming 
natural weighting in order to increase the Signal-to-Noise ratio in continuum and 12 CO. The 
resulting beamsize is of 0'.'61 x0'.'43. Only in the 2007 data C 18 showed detection above 3a 
and neither dataset detected 13 CO. 
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3. Results 
3.1. Continuum 

Compact 870/xm data shows two continuum peaks (Figured]). The unresolved stronger 
peak coincides with LMWOO's VLA1623A & B. A weaker peak located 10"5 to the west, 
matches with BA97's knot B. We will refer to this west peak with YC97's notation 
VLA1623W. Hereafter we refer to VLA1623A, VLA1623B and VLA1623W as A, B and W, 
respectively. 

Extended 870/im observations resolve A & B, but W is too weak to be detected in such 
short observations. B appears brighter than A possibly due to B being more compact than 
A and that much of the emission from A, present in the compact data, has been resolved 
out (Figure [Tb). 

The extended 1.3mm continuum observations resolved three peaks coinciding with the 
positions of A, B and W. A presents an elongated flattened structure perpendicular to the 
outflow axis, while B is compact and W is much dimmer (Figure WX)- 

We identify the infrared counterparts by comparing 870/xm and 1.3mm continuum 
to Spitzer data (Figure [Tk,b). Unresolved A & B are only detected in the MIPS1 band 
(40.4±4.24mJy) with the emission centered at A. Since B's SED starts to drop at 870/im 
(Figure |2]) we assume the MIPS1 emission is from A. Lack of detection of A & B in all 
IRAC bands indicates that both components are still deeply embedded and very young. 
Surprisingly we find that W coincides with a bright infrared source detected in the four 
IRAC bands and MIPS1 band. Resolved fluxes and uncertainties for each component are 
listed in Table [TJ 

We construct each component's SED using all available resolved fluxes from BA97, 
LMWOO, WKGAll and our work (Table [1]). The source size in all literature are in the 1~2" 
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range except in MIPS1 band; therefore, the flux we used originates from similar spatial 
scales. The MIPS1 detections are point sources, so their large source size is due to poor 
resolution. Fluxes between 24/im<A<7mm were fit by a single-temperature greybody fit 
with optical depth rcxz^ 3 . The best fit curves are shown in Figure [2] and fitting results are 
listed in Table [U The uncertainty in the absolute flux calibration may affect our fitting 



results. 



Ho et al 



f we adopt the most conservative estimate for SMA flux uncertainty (~30%. 



2004 ). the derived parameters (Tabled]) will only be altered by <10%. 



Using the dust temperature, To, obtained from the SEP fitti ng, assuming 



k 23 o = 0.01 cm 2 g 1 and a distance of 120±4.5pc (ILoinard et al. 



20081 ) . we derive the total 



(gas and dust) circumstellar mass M cir , column and number density (N# 2 & n H2 ) for each 
component from 1.3mm continuum (Table [TJ uncertainties were obtained through error 
propagation accounting for errors in distance, flux and Tp). Using the fluxes listed in 
Table[T] we derive the bolometric and submillimeter luminositites (L boi & L sufmm ) and their 
ratio for each component. We find L submm /L bo i = 1.17% a nd 0.18% for A and W, ratios 



chara cteristic of Class and Class I sources, respectively ( lYoung fc Evans 



2005 



Froebrich 



20051 ). For B we find L boZ = 8.9x10 5 L Q , considerably lower than VeLLOs. 



Dunham et al. 



( 120081 ) showed that incomplete SEDs lead to underestimated derived 
luminosities. To examine whether this affects A & W's classification, we calculate the 
luminosities from the best fit SED curve including the area under the extrapolated curve, 
L>sed, and then compare with the above derived luminosities (Tabled]). We find that L S ed 
are higher than the derived luminosities, as predicted, however the ratios are consistent 



with the evolutionary classification previously obtained. 



-8 - 



3.2. Outflow 



In Figure |3l we show VLA1623's outflows in 12 C0 (3-2) and (2-1) with a velocity range 
of -6-2 km s -1 and 4-15 km s -1 for the blueshifted and redshifted emission, respectively (v; sr 



=3.8kms 1 , YC97). Ove rlapped red- and 



transitions, as expected (lAndre et al. 



1990; 



ueshifted morph ologies are observed in both 



Dent et al. 



1995 



YC97), however we observe 



several additional features not previously reported. In 12 CO (3-2) the SE lobe is collimated 
with overlapped emission throughout, whereas the NW lobe presents a forked structure with 
overlapped red- and blueshifted emission only in the southernmost branch. In 12 CO (2-1) 
a slow blueshifted feature is located at W while a dim redshifted clumpy ouflow coincides 
with the northernmost branch of the NW lobe. In both transitions we observe a compact 
blueshifted knot elongated in the NE to SW direction and centered on B. The blueshifted 
knot presents a spatially overlapped redshifted counterpart only in 12 CO (2-1) (Figure [3]). 

SO spatially overlapped blueshifted (-0.5-2.3 km s -1 ) and redshifted (5-7.2 km s -1 ) 
emission as well as a systemic velocity component (3-4.4 km s" 1 ) and a southward 
blueshifted ridge are detected around A & B (Figure |3]). The SO emission is probably 
tracing circumstellar material shocked by A & B's outflows. 



3.3. Circumstellar Envelope 

Observed C 17 and C 18 emission centered on A present an extended flattened 
structure perpendicular to the outflow axis (Figure HJ), similar to A's shape in continuum. 
Both emission show similar velocity gradients characteristic of rotation along the major 
axis, seen in both velocity gradient maps and Position- Velocity (PV) diagrams (Figure H]). 
No velocity gradient is observed along the outflow axis, suggesting contamination from 
outflows is neglible. 
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C 17 0's dimensions (FWHM 370AUx240AU) are somewhat larger than A's 1.3mm 
continuum dimensions (240AUxl20AU). Assuming centrally dominated motion, we overlay 
Keplerian rotation curves (VocR -0 ' 5 ) for central masses of 0.1 and O.2M on the C 17 and 
C 18 PV diagrams. For comparison we overlay an infall curve VocR -1 . Although the infall 
curve passes through the peak of C 18 0, both emission are better traced by the Keplerian 
rotation curves. This suggests that both molecules trace the same rotating disk-like 
structure, with C 18 tracing the innermost region and C 17 the outermost region of the 
structure due to d ifferent resolutions. Th ough smaller, the characteristics of A's Keplerian 
disk are similar to 



Takakuwa et al. 



( 120121 ) 's Class I Keplerian disk, hinting at the beginnings 



of disk evolution. Neither emission appears associated with B and no circumstellar material 
was detected towards W. 



4. Discussion 

4.1. VLA1623A: The prototypical Class source 

From our SED fit of A we obtained T^ = 36K, higher than previously reported 
(Td = 15-20K, AWB93). This discrepancy may be due to our fitting using only resolved 
fluxes tracing the inner hotter regions of the core instead of considering the surrounding 
colder regions, as was done in AWB93. Despite the higher dust temperature, we find A 
to be a Class object based on lack of detection beyond 24//m and L submm /L bol >0.5%, 
in agreement with previous work. We do not consider the circumstellar-to-protostar-mass 
ratio as a criterion for evolutionary stage since we cannot account for all of A's envelope 
mass from our observations. 

A presents a flattened elongated envelope structure perpendicular to the outflow 
direction traced in C 17 and C 18 (Figure H]) as well as in continuum, most noticeable at 
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1.3mm (Figured]). Both C 17 and C 18 present signatures of rotation that suggest pure 
Keplerian rotation (Figure H]) and provide an estimate for the central object mass of A to 
be O.1-O.2M . Both the envelope traced by C 17 and C 18 and the 1.3mm continuum 
emission may be tracing a rotating disk about A, a possibility also raised by WKGAll. 
The presence of this rotating disk about A would be one of the first evidence of Keplerian 
disks around Class protostars, however further analysis must be carried out to confirm 
whether this is indeed a Keplerian disk. 



4.2. VLA1623W: A newly identified protostar 



We find W to be a true protostar, based on detection in all IRAC and MIPS1 
bands, as well as in the (sub)millimeter and centimeter regimes. We are unable to 
analyze Ws environment due to lack of detection of circumstellar material. Given the 
lack of circumstellar material emission towards W and the small circumstellar mass M cir 
derived from the 1.3mm continuum (Table [1]), we assume that Ws M cir is less than the 
central object mass M*. Additionally, AWB93 showed L su b/Lb o i<0.5% to be equivalent to 
M em ,<M,,,. Hence, based on an infrared spectral index am — 1.68±0. 07 (c2d catalogue), 
L-mhmm /Lhnz = 0.18%, M PM <M t and the Class-I-shaped SED, we classify W as a Class I 



source. 



Enoch et al. 



( 120091 ) mistook W as A+B due to their source identification criteria; 
however, from the listed source position and am the source is W, even more so since A 
lacks infrared e mission at A<24//m . The fact that W is a protostellar source confirms 



the suspicions of 



Dent et al. 



( 119951 ) and YC97 that the outflow and jet misalignments are 



caused by a second source. 
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4.3. The Nature of VLA1623B 



We find B to be an extremely cold (To ~5K) object with extremely low luminosity 
(~1O~ 5 L ). No associated C 17 and C 18 emission suggest strong CO depletion, consistent 
with the derived extremely cold dust temperature (CO sublimation temperature ~20K). 
We further suspect that B may be driving a pole-on outflow (see section I4~4l and Figure [3]). 
While we cannot determine B's nature unambiguously, some possibilities can be ruled out. 



Caratti o Garatti et al 



( 120061 ) show VLAl623's jet knots have gas temperatures of 
1700-4200K and velocities between 60-100 km s -1 . The strong CO depletion at B indicates 
that the gas temperature is <20K, far below that e xpected for jet knots. Observations 



20121 ) show obvious 



of monopolar continuum jet knots in T Tauri stars ( [Rodriguez et al. 
shifts in the knot positions in short periods of time. Assuming simple ballistic motion and 
velocities of 60-100 km s -1 , B should have shifted from BA97's initially reported position by 
~2-3" in a period of 15 years. However, comparison of BA97's and our 2012 observations 
show no such shift. VLA high re solution (~0f!5) centi meter continuum observations of 



jet-driving Class and I objects (IReipurth et al 



20041 ) found most detected continuum 



peaks to be protostellar sources with continuum jets seen as extended emission in the 
direction of the source's outflow/jet. B is observed as a separate compact component and 
is coincident with counterparts in the (sub)millimeter regimes, making it more like a YSO 
and less like a jet knot. Although VLAl623's alignment and distribution might resemble 
a jet, several MPSs are observed to have th eir components e ither align ed or within their 



companion's outflow (e.g. NGC1333 SVS13: 



Chen et al 



2009 



L1448 C: 



Hirano et al. 



20101) 



possibly caused by filament fragmentation. Thus VLA1623's configuration is not highly 
unlikely. Furthermore, a quick Jeans Mass calculation shows that given B's mass, dust 
temperature and density, the core is unstable and most likely collapsing. Hence even if it 
were product of outflow-environment interaction it is on its way to becoming a YSO. Thus, 
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B cannot be a jet/outflow feature as proposed by MOA12. 



Other possibilities for the nature of B are a starless core or proto-brown dwarf. While 



the CO deplet ion and extremely 



starless cores ( Tafalla et al. 



2002 



ow temperature an d luminosity are characteristics of 



Aikawa et al. 



20051 ). the mass, density and presence of 



outflow rule out the possibility of B being a starless core. If we assume a star formation 
efficiency of 30% the resulting object from B would be a ~0.1M Q star, well above the brown 
dwarf limit. With current data and B's physical and chemical characteristics, we cannot 
determine the eventual outcome but we can conclude that B is an extremely young object 
between the starless core and Class stages. 



4.4. A triple system 

Comparison of both 12 CO maps and the presence of three sources suggest that 
VLAl623's outflow is not product of a single source. We propose that the NW lobe is not 
a cavity but instead two outflows. A drives the northernmost branch while W drives the 
southermost branch based on its position and the 12 CO (2-1) blueshifted feature. The SE 
lobe is difficult to disentangle from the observations presented here, requiring comparison 
of 12 CO and jet observations. 

For B, the observed spatially overlapped 12 CO (2-1) emission suggests a compact 
pole-on outflow. The red- and blueshifted spatially overlapped SO emission and absence of 
a redshifted lobe in 12 CO (3-2) further support this morphology, since the redshifted lobe of 
B would be located directly behind the blueshifted lobe which is stronger and more easily 
detected by the compact configuration observations. Additionally, the 12 CO (3-2) redshifted 
lobe of B is easily lost among the much stronger redshifted emission of the outflows from A 
& W. 
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Further analysis on VLAl623's kinematics, outflows and whether it is gravitationally 
bound, as well as determining B's nature, are required to fully disentangle VLA1623. 
However, from its alignment, physical characteristics, environment and differing evolutionary 
stages we can conclude that VLA1623 is a triple non-coeval system probably formed 
through the delayed fragmentation of a filament. These characteristics make VLA1623 a 
good target for constraining and studying fragmentation in low-mass protostars. 
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Fig. 1. — Continuum: Spitzer (a) IRAC4 (8/xm) and (b) MIPS1 (24/um) halftone images; (c) 
Extended 870 /mi continuum; (d) Very Extended 1.3mm continuum with alternate component 
names. All images are overlaid with 870/mi compact continuum contours (light grey in (c) 
and (d)) for scale comparison. Dashed and dashed-dot circles show the SMA primary beam 
for 870/zm and 1.3mm, respectively. RMS noise (mJybeam -1 ) is in curly brackets. Contours 
are in steps of -5a, -3cr, 3cr, 5a, 10a, 30a and 60a for 870/rni compact continuum; -4a, -3a, 
3a, 4a and 5a for 870/im extended continuum; and -5a, -3a, 3a, 5a, 10a, 20a and 50a for 
1.3mm Very extended continuum. (Color figure in the online journal) 
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Fig. 2. — SED for each component of VLA1623 constructed using data from BA97, LMWOO, 
WKGAll and this work (Tabled]). The solid black line shows the best single-temperature 
graybody fit for fluxes at 7mm>A>24/mi, with the best fit dust temperature and dust 
emissivity index (3 written in the upper right corner. A simple interpolation (dashed line) is 
made for points outside the fitting range. 
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Fig. 3.— Intensity integrated maps for SO (9-8), 12 C0 (3-2) and 12 C0 (2-1). Redshifted 
and blueshifted emission are shown in black and grey contours, respectively. RMS noise 
(Jybeam -1 kms -1 ), integrated blue- and redshifted velocity ranges (kms -1 ) are listed in 
curly, round and square brakets, respectively. Contours are in steps of 10a, 20cr, 30a, 40a, 
50a, 90a, 130a and 170a for 12 CO (3-2); 4a, 6a, 7a, 12a and 25a for 12 CO (2-1); and 3a, 6a, 
12a and 18a for SO (9-8). Negative contours not plotted for clarity. Positions of VLA1623A, 
B & W are marked with a square, triangle and circle, respectively. (Color figure in the online 
journal) 
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Fig. 4. — Intensity integrated (contours) and Velocity gradient (halftone) maps for (a) 
C 17 and (b) C 18 0. The diagonal line in both maps indicates the Position- Velocity (PV) 
diagram cut (P.A. = 30°) for (c) C 17 and (d) C 18 0. RMS noise (Intensity integrated: 
Jybeam -1 kms -1 ; PV diagram: Jybeam -1 ) is listed in curly brakets. Intensity integrated 
map contours are in steps of -4a, 4a, 7a, 11a and 13a for C 17 0; and -2a, 2a, 3a, 4a, 5a, 
6a and 7a for C 18 0. PV diagram contours are in steps of -3a, 3a, 5a, 7a, 9a, 11a and 13a 
for C 17 0; and -2a, 2a, 3a, 4a and 5a for C 18 0. Keplerian rotation curves (VocR -0,5 ) for 0.1 
and 0.2M© and infall curve (VocR -1 ) are plotted on the PV diagrams. The vertical and 
horizontal solid lines in the PV diagrams show v^ r = 3.8kms _1 and VLAl623A's position, 
respectively. VLA1623A and VLA1623B are marked by a square and triangle, respectively. 
(Color figure in the online journal) 
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Table 1. Observations and Physical Parameters of VLA1623 





VLA1623A 


VLA1623B 


VLA1623W 


Ref. 


RA 


16:26:26.397 


16:26:26.308 


16:26:25.628 


1 


Dec 


-24:24:30.71 


-24:24:30.57 


-24:24:29.60 


1 


Separation (") a 





1.18 


10.51 


1 


Integrated Flux densities (mJy) b 








1.9 ± 0.09 


2 


4.5/^m 






6.6 ± 0.2 


2 


5.8//m 






0.5 ± 0.7 


2 


8/xm 






7.2 ± 0.2 


2 


24yum 


40.4 ± 4.2 C 




272.2 ± 15.1 


2 


870/rni 


1153 ± 100 d 


87.2 ± 12 


128 ± 16 


1 


1.3mm 


152 ± 8 


106 ± 8 


25.7 ± 8 


1 


2.7mm 


25.5 ± 6.3 


25.8 ± 3.5 


3.1 ± 1.5 C 


3 


7mm 


1.2 ± 0.2 


2.0 ± 0.2 




4 


3.6cm 


0.1 ± 0.02 


0.4 ± 0.02 


0.2 ± 0.02 


5 


4.4cm 


< 0.2 


< 0.2 




4 


6cm 


< 0.4 ± 0.02 


< 0.4 ± 0.02 


0.2±0.02 


4 



SED fitting Results 
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Table 1 — Continued 



VLA1623A VLA1623B VLA1623W Ref. 



Td (K) 


36 ± 5 


5 ± 1 


41 ± 5 


1 




1.2 ± 0.6 


1.7 ± 0.5 


1.6 ± 0.4 


1 


Reduced % 2 


21.7 


5.8 


0.7 


1 


Physical Parameters 


Lfeo/ (L©) 


0.03 ± 0.003 


(8.9 ± 0.5) x 10~ 5 


0.03 ± 0.002 


1 


^->submm (L©) 


(3.9 ± 0.3) x 10" 4 




(4.8 ± 0.5) x 10~ 5 


1 


Ls«6mm/Lfeo« (%) 


1.2 ± 0.1 




0.2 ± 0.02 


1 


L>bol,SED (L©) 


1.1 ± 0.2 


(2.9 ± 0.1) x 10" 4 


1.4 ± 0.2 


1 


Ls«6mm,SEO (L©) 


(1 ± 0.6)x 10" 2 




(5 ± 2) x 10- 3 


1 


Ls«6mm,S'BD/L6oi,S , £;D (%) 


0.9 ± 0.5 




0.4 ± 0.2 


1 


M cir (M©) f 


0.02 ± 0.004 


0.3 ± 0.2 


0.003 ± 0.001 


1 


N H2 (cm- 2 ) f 


(4.4 ± 0.8) x 10 23 


(9.9 ± 5.0) x 10 24 


(3.6 ± 1.2) x 10 23 


1 


n# 2 (cm~ 3 ) f 


(2.2 ± 0.4) x 10 8 


(5.5 ± 2.8) x 10 9 


(4.1 ± 1.4) x 10 s 


1 



Reparation with respect to VLA1623A 



b Flux uncertainties are statistical 

c Unresolved VLA1623A & B flux. However, since VLA1623B peaks at around 1000/im, 
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its contribution to the 24/im flux is neglible. 

d Flux calculated by subtracting VLA1623B's 870/um extended integrated flux from the 
unresolved VLA1623A & B 870/im 870/im compact flux. 



e We obtained reduced BIMA data from 
selves. 



Looney et al. 



(120001 ) and calculated the flux our- 



E Calculated from 1.3mm continuum emission 



References. — (l)This work; (2)c2d Catalogue; (3)LMW00; (4)WKGA11; (5)BA97 



